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Abstract 



By variational Monte-Carlo method developed Ceperley et al. for the sim- 
ulation of fermi systems in macroscopic conhning potential well we simulate 
various spin ground states of the coulomb clusters with 2,3 and 4 particles 
in a broad two-dimensional {2D) parabolic well. In this method quantum 
state numbers determining the variational wave function are not the numbers 
of well quantum states but numbers of the equilibrium spatial positions of 
particles that give a minimum of the system potential energy. The ground 
states with parallel, antiparallel spins and as well, as bose state are simulated. 
For the cluster with three particles it is also simulated the state when two 
particles have one direction of spin and third opposite. The simulation shows 
that clusters with parallel spins have lower ground state energy than clusters 
with other spin configurations and bose state. That reminds a Hund's rule in 
atomic physics when in not full filled atomic shells electrons prefer to have a 
state with parallel spins. 



I. INTRODUCTION 

A considerable experimental succes in localization and cooling of several ions in electro- 
static traps and recent technological advances in the fabrication of the quantum dots in 
semiconductors have induced a big theoretical interest [Ij-^] to these objects. A 2D 
parabolic potential well is a good model potential well for the confining of electrons in the 
quantum dots ||^ in many cases. In the works it was shown the equivalence of coulomb 
system in effective electrostatic potential Uext = o,R^ to the Tomson model of atom. The 
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classic and quantum properties of coulomb clusters in traps were also investigated in work 
i- 

Here we consider the coulomb clusters that correspond to quasi two dimensional systems 
of ions or electrons taking place in strong anisotropic 2D traps that can be supposed as 
broad. By variational Monte-Carlo method developed Ceperley, Chester and Kalos |Ty] for 
the simulation of fermi systems in macroscopic confining potential well we simulate various 
spin ground states of the coulomb clusters with 2,3 and 4 particles in a broad 2D parabolic 
well. In this method quantum state numbers determining the variational wave function 
are not the numbers of well quantum states but numbers of equilibrium spatial positions 
of particles that give a minimum of the system potential energy. The ground states with 
parallel, antiparallel spins and as well, as bose state are simulated. For the cluster with 
three particles it is also simulated the state when two particles have one direction of spin 
and third opposite. The simulation shows that clusters with parallel spins have lower ground 
state energy than clusters with other spin configurations and bose state. That reminds a 
Hund's rule in atomic physics when in not full filled atomic shells electrons prefer to have a 
state with parallel spins. 



II. MODEL OF SYSTEM AND THE SIMULATION METHOD 

Let us consider a system - cluster having particles in a 2D parabolic well with Hamil- 
tonian 

^ = -A^ + f/po*(i?), (1) 

where R corresponds to the two dimensional vectors (ri, r2, tat) and the potential energy 
is 

N N 

U,ot{R) = Y.^.^ + Y.\^^-r,\~' (2) 

i=l i>j 

Hamiltonian (|T]) contains a single quantum parameter A = fi^ a^^^ (2me^^^)~^ , where m 
and e mass and charge of coulomb particle respectively, parameter a determines a parabolic 
well Uext = o.R^ which locates the particles. For a broad 2D parabolic well we assume A < 1. 



We simulate the coulomb clusters by variational Monte-Carlo method proposed in |T0 
for the simulation of fermi systems in macroscopic confining potential well. It is based on 
the trial variational wave function 

^{R) = MR)D, (3) 

where to take into account an interparticle electrostatic correlation and a 2D parabolic well 
it is introduced the Jastrow wave function: 

MR) = exp ( -/3Er-;2 - 7E l^"'^ - ^".rM • (4) 
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The expression for the quantity D depends from the system simulating state. If it is simu- 
lating a bose state of the system then 



N 

D = nexp[-a(f,-r;o)']. (5) 

i=l 



If there is a state of pure fermions then in according with ||T0 

a 



D=X{dei{Dt^), (6) 



=1 



where spin index s has g number of states (so for the parallel spin state g = 1 and for the 
antiparallel spin state g = 2) and det{D^j) is Slater determinant, where 

Dt^ = exp[-aif, - f,o)%. (7) 

In the expressions @ and (j^) the vectors fio provide a minimum of the potential energy 
Upot{R) (H), i-e. a minimum of the potential energy for the classical system. In the expression 
(0) index i indicates the particle's coordinate and j indicates a quantum state. 

The antiparallel spin state for the cluster with two particles coincides with bose state of 
system. 

The antiparallel spin state for the cluster with four particles has a product of two Slater 
determinants. We assume that in (|^ j = 1 and j = 4 describe an one direction of spin and 
j = 2 and j = 3 the opposite direction. Therefore, in the expression first determinant 
determines by j = 1,4 and second by j = 2, 3. 

The paralell spin state for the cluster with four particles has only one Slater determinant, 
where j = 1, 2, 3, 4. 

A minimum of the ground state energy 

^ !dm*iR)H^iR) 

IdR\^{RW ^ ^ 

is achieved by variation of parameters a, (3 and 7 in the expressions (il)-®. It is assumed 
that variational parameters (3 and 7 in the Jastrow wave function (^) are the same for all 
spin and bose states of all clusters (see |jlO|)- That means that there is no influence of the 
interparticle electrostatic correlation and external field on the spin part of the wave function. 

The cluster ground state simulation is occured by Metropolis algorithm. In this algorithm 
each particle of the cluster is tested on the uniform random displacement inside square with 
side A from center of square VjQ. New spatial position of the particle is accepted by the 
probability 

P = min[l,|^(r„eJ/V'(r)r]. (9) 

The size of A is chosen by the such way that the ratio of accepted number to whole number 
of particle displacements has been roundly equal 1/2 for all particles. Thus, the full variation 
of energy (P) is taken place by the variation of four parameters a, /5, 7 and A. 
The mean values of the quantities in algorithm are calculated by the formula: 
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J dRiP*{R)F{R)iP{R) 



1 M 



(10) 



jdRmRw 

where Ri is the i-th spatial configuration of the system whole number of which is equal M. 

The full energy of the system is the sum of every particle's energy. Let's consider the 
expressions for the kinetic and potential energies, for example, for the first particle. 

The kinetic energy is 

^*^[(Vi In + In ^]dR. 



2m 



Erin - 

If we introduce the determination 

Ti 



(11) 



V?ln^ 

i 4m 



(12) 



then the expression for the full energy will have a form: 



where 



N g2 



ar 



i=2 



(13) 
(14) 



In the expressions (|T2|) and ( p^ the angle brackets mean the mean value (p^Ol). 

In the paper it is shown that Ti = Fl and therefore, we are calculating the energy 
El by the formula 

Ei = Vi + Fl (15) 
For the bosons containing cluster the expression for F^ is 



^ 2m 



N 



- E Viu(rij) - Vix(ri) + Vi$i(ri) 

Here u{rij) = 'y/{rj — ri), x(ri) = Prf and 

$i(ri) = exp[-a(fi - fio)^]. 
If it is simulated the fermions containing cluster then 



(16) 



^ 2m 



N N 

Y: ^Mn,) - Va(ri) + E D^i^i^An) 

j=2 j=l 



(17) 



Here D,- is the inverse matrix of matrix Df,. 

The main complexity of the simulation by the variational Monte-Carlo method of big 
number fermions containing systems is connected with technical difficulty of inversion of 
matrix D^y In the paper |T^] is proposed a simple algorithm for the inversion of matrix Df^- 
for these systems. Our clusters have a small number of particles. Therefore, we use direct 
method of determination of elements of inverse matrix by the elements of direct matrix 
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III. THE SIMULATION RESULTS 



As it was said above, Hamiltonian (0) has a single quantum parameter A = 
h'^ a^^^ {2me^^^)~^ . In accordance with this determination of parameter A all length quantities 
are expressed in the e^^^/a^/^ units and energy quantities in the e'^^^a^^^ units. 

The positions of particles ro that give a minimum of the potential energy Upot{R) (@) ^^i 
the classic systems of the clusters are determined analytically. For all spin and bose states 
of one cluster they are suggested identical. 

For the cluster with two particles the positions of tq situate on the one line, on the 
distance 1/2 from the centre of coordinates that coincides with the centre of a parabolic 
well. 

For the cluster with the three particles the coordinates of ro are the coordinates of the 
triangle 

r - (^V^'^- r -r^V^'l. r - Q- 
roxi — — \^j 2' ^""^ ~ V 2/ 2' ~ ' 

3\V3 1 /3xi/3 1 /3\V3 1 



For the cluster with four particles a minimum of Upot{R) provides by square (with par- 
ticles on the vertices of square and the centre of square on the centre of a parabolic well 
and coordinates). The distance between the centre coordinates and square vertex is equal 

ro=((4 + v^)/(8v^))'^'. 

Variational parameters (3 and 7 that are the same for all spin and bose states of all 
clusters we determine by the simulation of the bose state for the two particle cluster. The 
numerical quantities of these are f3 = 0.5 and 7 = 1.5. 

Whole number of displacements inside A for every particle we choose 10"^, so absolute 
exactness of the calculated mean energies is roughly equal 10~^ ( we would like to note that 
for the Markov random processes, i.e. for independent random processes, absolute exactness 
of the calculated in simulation mean values is proportional to 1/ \/N, where is number of 
random simulations). 

There is no a necessity for the simulation of the energies at every numerical quantity of A 
for every spin or bose state of the one cluster. Let we have the calculated in the simulation 
the potential energy V and the kinetic energy Ekin at any one fixed A. As the kinetic energy 
in Hamiltonian (|I]) is directly proportional to A, so for the getting of the kinetic energy for 
another Ai one can recalculate it by the formula Ekin{.^i) = E^inW ■ Ai/A; the potential 
energy V does not depend from A and one can be calculated once for this spin or bose state 
of cluster. 

Fig. 1 shows variational parameter a dependence (in the e^^^a^/^ units) of the potential 
energies V^bs, V^afs and V^ps ( bose, antiparallel spin (antiferrospin) and parallel spin 
(ferrospin) states respectively) for the cluster with four particles, of the full energies E2BS 
and E2FS ( bose and ferrospin states respectively) for the cluster with two particles and of the 
kinetic energies Ekinabs^ Ekinaafs and Ekinafs per one particle at the fixed parameter 
A = 0.01. 

In this figure the potential energy for the bose state V^ss is increased when parameter 
a goes to zero. This takes place due the increasing of the electrostatic interaction between 
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coulomb particles because the decreasing of a means the increasing of the mean size of the 
wave function per one particle. The a dependence of the potential energy for parallel spin 
state V4FS when a goes to zero has opposite character. This is occured due the increasing 
of the exchange interaction between fermions. 

Fig. 2 shows the A dependence (in the e^^^a^/^ units) of full ground state energy for every 
cluster per one particle. In this figure the symbol 3AFS corresponds to the state when two 
particles have one direction of spin and third opposite. The numerical quantities of the full 
energies in Fig. 2 are the minimums of these energies on variational parameter a. 

At last, for the demonstration of qualitative behaviour of the module of the wave function, 
we have simulated it for the cluster with three particles at parameters P — 1.7 and 7 = 1.0 
(though these parameters of f3 and 7 do not give the minimum of the full energy, but they 
provide for the module of the wave function an obvious visual illustration). This quantity 
is outlined in the figures 3,4 and 5 for the bose, parallel spin and mixed spin ( i.e. 3AFS) 
states respectively. From these figures one can see clearly that exchange interaction between 
fermions moves apart a quantum spatial distribution of the cluster particles. Bose and fermi 
like together behaviours of the module of the wave function are outlined in Fig. 5. 

The main result of the present work is outlined in the figures 1 and 2. It turns out that for 
the every cluster with 2,3 and 4 coulomb particles in a 2D parabolic well the ground states 
with parallel spins are energetically advantageous than other spin and bose states. That 
reminds a Hund's rule in atomic physics when by taking account of exchange interaction in 
not full filled atomic shell it is profitable for the electrons to have a single direction of spin. 

We are grateful Yu.E. Lozovik for indication on this task and helpful discussions. 
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FIG. 1. The a dependence of the potential energies V4BS, Vaafs, Vafs, the full energies E2BS, 
E2FS a^nd the kinetic energies Ekinabs^ Ekinaafs^ Ekinafs at the fixed A = 0.01. 
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FIG. 2. The A dependence of the full energies for all clusters. 



FIG. 3. The module of the wave function for the bose state for cluster with three coulomb 
particles in a 2D parabolic well. 




FIG. 4. The module of the wave function for the parallel spin state for cluster with three 
coulomb particles in a 2D parabolic well. 
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FIG. 5. The module of the wave function for the mixed spin state for cluster with three 
coulomb particles in a 2D parabolic well when two particles have one direction of spin and third 
opposite. In the figures 3,4 and 5 the numerical quantities in the vertical axes expressed in single 
numerical units. 
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